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The optical absorption of the Frohhch polaron model is obtained by an approximation-free Di- 
agrammatic Monte Carlo method and compared with two new approximate approaches that treat 
lattice relaxation effects in different ways. We show that: i) a strong coupling expansion, based on 
the the Franck-Condon principle, well describes the optical conductivity for large coupling strengths 
(a > 10); ii) a Memory Function Formalism with phonon broadened levels reproduces the optical 
response for weak coupling strengths (a < 6) taking the dynamic lattice relaxation into account. In 
the coupling regime 6 < a < 10 the optical conductivity is a rapidly changing superposition of both 
Franck-Condon and dynamic contributions. 
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The Franck-Condon (FC) principle(lj] and its validity 
have been widely discussed in studies of optical transi- 
tions in atoms, molecules^ and solidsQ- For instance, 
studies on the dynamics governing color center spec- 
tra, Mossbauer spectra^, and tunneling in polarizable 
mcdiaQ have established the limits of applicability of 
the FC principle in solids. Generally, the FC princi- 
ple means that if only one of two coupled subsystems, 
e.g. an electronic subsystem, is affected by an exter- 
nal perturbation, the second subsystem, e.g. the lattice, 
is not fast enough to follow the reconstruction of the 
electronic configuration. In the opposite limit, when the 
perturbation is slow or even static, the characteristic time 
of lattice interconfigurational coupling Tic is short enough 
for the lattice to follow slow changes of the electronic 
state dynamically. It is clear that the justification for 
the FC principle is the short characteristic time of the 
measurement process T^p ^ Tic, where Tmp is related 
to the energy gap between the initial and final states, 
AE, through the uncertainty principle: T„ip — h/{AE). 
Then, the spectroscopic response considerably depends 
on the value of the ratio Tmp/ Tic- For example, in mixed 
valence systems, where the ionic valence fluctuates be- 
tween the configurations and with characteristic 
time Tic ~ 10~^'^s, the spectra of fast and slow experi- 
ments are dramatically different Photoemission ex- 
periments with short characteristic times Tmp ~ 10~^^s 
(FC regime), reveal two lines, corresponding to and 
states. On the other hand slow Mossbauer isomer shift 
measurements with T^p ~ 10~^s show a single broad 
peak with mean frequency between signals from pure 
P and shells. Finally, according to the paradigm 
of measurement process time, magnetic neutron scatter- 
ing with T„ip ~ Tic revealed both coherent lines, with 
all subsystems dynamically adjusted, and broad incoher- 



ent remnants of strongly damped excitation of and 
shells |3- Actually, the meaning and the definition of 
the times Tic and T^p varies with the system and with 
the measurement process although the spectroscopic re- 
sponse essentially depends on whether one of the two in- 
teracting subsystems is adjusted to the changes, induced 
by a measurement probe, in the other subsystem. 

To investigate the interplay between the measurement 
process time Tmp and the adjustment time t,;c, we study 
in this Letter the optical conductivity (OC) of a paradig- 
matic model for electron-phonon (e-ph) interaction: the 
Frohlich model. Our aim is to investigate the OC from 
the weak to the strong coupling regime (in this model 
the e-ph coupling strength is controlled by the dimcn- 
sioless parameter a) by three methods: i) the diagram- 
matic Monte Carlo (DMC) method |^ 0| which gives 
numerically exact answers in all e-ph coupling regimes; 
ii) the memory function formalism (MFF) which is able 
to take dynamical lattice relaxation into account; iii) and 
a strong coupling expansion (SCE) which assumes the 
FC principle. In the coupling regime 6 < a < 10 we find 
that the OC spectrum exhibits two features with differ- 
ent behavior. The higher-frequency feature quickly de- 
creases its spectral weight with increasing coupling con- 
stant whereas the lower-frequency feature does the op- 
posite. Besides, the numerically exact calculations of the 
OC (DMC) follow the prediction of the extended MFF 
for a < 6, while they are in fair agreement with SCE 
for a > 10. We conclude that nonadiabaticity destroys 
the FC classification for a < 10 while the FC principle 
rapidly regains its validity at large coupling strengths due 
to the fast growth of the energy separation between the 
initial and final states of the optical transitions. Fur- 
thermore, both adiabatic FC and nonadiabatic dynami- 
cal excitations coexist in the intermediate e-ph coupling 
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regime, 6 < a < 10. The crossover is controlled by the 
adjustment time Tic ~ ^|^^^ set by the typical nonadia- 
batic energy V (see text below). 

In the Frohlich polaron model^Jl the electron (r and 
p are the position and momentum operators) is scat- 
tered by phonons (at denotes the creation operator with 
wave number q) with c-ph interaction vertex Mq ~ 
ihuo (i?p47ra/g2y)^/^ 

H = p^/2m + hujQ ata,- + ^[M,e*«-%- + h.c]. (1) 
Here a is the dimensionless coupling constant, Rp = 

1/2 

{h/2muJo) , and V is the volume of the system. The 
band mass of the electron to, Planck's constant ?i, the 
electron charge and the dispersionlcss longitudinal opti- 
cal phonon frequency ujq are set below to unity. Although 
the OC of this model attracted attention for vears|l2j. 
there exists no analytic approach giving a satisfactory 
description for all coupling regimes. The most success- 
ful approach is that based on the Feynman path inte- 
gral technique [13 (DSG), where the OC is calculated 
starting from the Feynman variational model (FVM) [l^ 
for the polaron and using the path integral response 
formalism^3- Subscquentlythe path integral approach 
was rewritten in terms of the memory function formal- 
ism (MFF)[il]. These approaches are completely suc- 
cessful at small e-ph couplings, are able to identify some 
of the excitations at intermediate and strong e-ph cou- 
plings without reproducing the broad structures present 
in DMC data 0. 

Extended memory function formalism. In order 
to solve the aforementioned problem regarding the 
description of the OC main peak line width at in- 
termediate e-ph couplings, we modified the DSG 
approach to include additional dissipation processes, 
whose strength is fixed by an exact sum rule. Within 
mff|i3 the interaction of the charge carriers with 
the free phonon oscillations can be expressed in terms 
of the electron density-density correlation function, 
x{Q,t) = —iOU) (exp \iq- r{t)]exv \—iQ- r(0)]), which is 
evaluated in a direct wayjl^ using the FVM, where 
the electron is coupled via a harmonic force to a 
fictitious particle that simulates the phonon degrees of 
freedom. Within this procedure the electron density- 
density correlation function takes the form: Xm('f, ^ 
-ie{t)exp [-iqH/2M] exp [-q^R{l - e-™*)/2M] , 
where R = (A/ — l)/t;, M (the total mass of electron and 
fictitious particle), and v are determined variationally 
within the path integral approach [l^. The associated 
spectral function Amiq^oj) = — 2Qxm(9, i^) is a series of 
5 functions centered at q^ jlM + nv (n is integer). Here 
q^ /2M represents the energy of the center of mass of 
electron and fictitious particle, and v is the energy gap 
between the levels of the relative motion. To include 
dissipation we introduce here a finite lifetime for the 



states of the relative motion, which can be considered as 
the result of the residual e-ph interaction not included 
into the FVM. To this end, in XmiQ^t) we replace the 
factor exp[—ivt] with (1 + it/T)~^^ which leads to 
the replacement of i5 functions by Gamma functions 
with mean value and variance given respectively by 
q^ /2M + nv and nv/r. The parameter of dissipation r 
is not an adjustable parameter but is determined by the 
third sum rule for A[q, w), which is additional to the first 
two ones that are already satisfied in the DSG model 
without damping. As expected, r turns out to be of the 
order of lOq"^ . If broadening of the oscillator levels is 
neglected, t ^ oo, the DSG resultsHHIi 

are recovered. 

Within this latter approach [l^ llij the polaron OC was 
previously interpreted in terms of relaxed excited states 
transitions, FC-transitions and transitions to scattering 
states |il|l3- 

Strong coupling expansion. In the limit of strong cou- 
pling strengths the adiabatic lattice deformation is large 
and the lattice kinetic energy can be regarded as a per- 
turbation. Within the Landau and PekarfjfijfLP') ap- 
proach adiabatic lattice displacements arc taken into ac- 
count by a unitary transformation and the Hamiltonian 
is divided into two contributions. H = Hq + Hj: the first 
one describes the electron oscillation in a self-consistent 
quadratic potential with frequency = Aa^/Qn and 
the second one represents the residual e-ph interaction. 
To get a quantitative estimate for the characteristic fre- 
quency of the quadratic potential, we improve the LP 
value considering the effects of the translational invari- 
ance and the residual interaction iJr |2 ij. It turns out 
that the frequency uji = (4Q;^/97r — 3.8) differs from that 
of LP Lu^^ by a constant shift 3.8. Then, starting from 
the Kubo formula, taking into account all multiphonon 
processes and neglecting recoil as well as correlation be- 
tween the emission and absorption of successive phonons 
in all orders of the perturbation Hj, one arrives to the 
following expression for the real part of the OC 

5Ro-(w) = o-Qt^ } — ;- i^^sT 5[uj - UJ2 - n], (2) 

resembling that expected for an exactly solvable indepen- 
dent oscillators model Parameters, specific for the 

Frohlich polaron model, are 0^2 = (j^i — Ws, (Tq = tt/2uji 
1/2 

and ujs ~ a {uJi/16n) . The parameter uji is the FC 
transition energy, and uJs is the energy shift due to lat- 
tice relaxation. Naturally, for large enough coupling 
strengths the envelope of the Poisson distribution (0) is 
well described by the Gaussian 

^'^(")="(W^^"P1 2^j- 

Interpretation of the DMC results. As expected, in the 
weak coupling regime (Fig. la), both our extended MFF 
with phonon broadening and DSG|l3j are in very good 
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FIG. 1: Comparison of the optical conductivity calculated within the DMC method (circles), extended MFF (solid line) and 
DSG0,01 (dotted line), for four different values of a. The arrow indicates the lower-frequency feature in the DMC data. 




U) to U) a 

FIG. 2: a), b) and c) Comparison of the optical conductivity calculated within the DMC method (circles), the extended 
MFF(solid line) and SCE (dashed line) for three different values of a. d) The energy of the lower- and higher- frequency 
features obtained by DMC (circles and triangles, respectively) compared with the FC transition energy as calculated with the 
SCE (dashed line) and with the energy of the peak obtained from the extended MFF (solid line). In the inset the weights of 
Franck-Condon and adiabatically connected transitions are shown as a function of a. We have used for -q the value 1.3. 



agreement with the DMC data|9|, showing significant im- 
provement with respect to the weak couphng perturba- 
tion approach which provides a good description of 
the OC spectra only for very small values of a|2J|. For 
4 < a < 8, where DSG underestimates the peak width 
(Figs, lb to Id), the damping, introduced in the extended 
MFF approach, becomes crucial. Results of the extended 
MFF are accurate for the peak energy and quite satisfac- 
tory for the peak width (Fig. Ib-d). 

However, we observe a very interesting and unexpected 
behavior in the coupling regime 6 < a < 10. Two fea- 
tures are present in the OC given by DMC: the position of 
the lower-frequency peak (or shoulder) corresponds to the 
predictions of the SCE (eas l2l3|l while that for the higher- 
frequency peak follows the extended MFF value (Fig. 2a). 
The higher-frequency feature rapidly decreases its inten- 
sity with increasing a and, at large values of a (Figs. 2b 



and 2c), the OC given by DMC is in fair agreement with 
the SCE results, which are strongly dominated by the 
FC transitions. Finally, comparing the peak and shoul- 
der energies, obtained by DMC, with the peak energies, 
given by MFF, and the FC transition energies from the 
SCE (Fig. 2d), we conclude that as a increases from 6 to 
10 the spectral weights rapidly switch from the dynamic 
regime, where the lattice follows the electron motion, to 
the adiabatic regime dominated by FC transitions, where 
the nuclei are frozen in their initial configuration. 

Breakdown of the FC picture. In order to support this 
scenario we present an analytical estimation of the FC 
breakdown based on the following arguments. The ap- 
proximate adiabatic states Xi.;9((5)V'i(r, Q), where i is 
the electronic index and Xi,l3{Q) is the eigenfunction of 
the lattice connected with the electron wave function 
ipi{r,Q), are not exact eigenstates of the system. These 
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states are mixed by nondiagonal matrix elements of the 
nonadiabatic operator V and the exact eigenstates are 
linear combinations of the adiabatic wave functions. Be- 
ing interested in the properties of transitions from the 
ground {G) state to an excited (EX) state, whose en- 
ergy corresponds to that of the OC peak, we consider 
mixing of only these states and express the exact wave 
functions '^g,ex{y,Q) as linear combinations [j^l 

^G,EX = ^g'/3 Xg,l3'4'g+^ex,!3' Xex,ii'^ex of the adiabatic 
groundstate {Xg,p4'g) and the adiabatic excited state 
(Xes./s'V'ea:)- The supcrposition coefficients are deter- 
mined from standard techniques |25ll26j where the nondi- 
agonal matrix elements of the nonadiabatic operator |25j 
are expressed in terms of matrix elements of the kinetic 
energy operator M, the energy gap between the excited 
and ground states Ai? = E^x — Eg, and the number np 
of phonons in the adiabatic state: 



whether the lattice can or cannot follow adiabatically the 
change of electronic state at the transition. 

Conclusions. Comparing numerically exact data on 
the optical conductivity, obtained by the Diagrammatic 
Monte Carlo method, with results from our Extended 
Memory Function approach with phonon broadened lev- 
els, which takes dynamic lattice relaxation into account, 
and results of the strong coupling expansion we found 
that the Franck-Condon picture breaks down at a < 10. 
The breakdown of the Franck-Condon picture is caused 
by nonadiabatic mixing of initial and final states, which 
destroys the Franck-Condon classification scheme and, 
hence, the excitation processes with dynamic adjustment 
of the lattice start to dominate. Finally we find evidence 
for an intermediate coupling regime 6 < a < 10 where 
static and dynamic lattice responses coexist. 



P± = M{AE)-\ Ufj + 1/2 ± 1/2 + M^iAE)- 



(4) 



The extent to which the lattice can follow a transition 
between electronic states, depends on the degree of mix- 
ing between initial and final exact eigenstates through 
the nonadiabatic interaction. If the states are strongly 
mixed, the adiabatic classification has no sense, the FC 
transitions have no place and the lattice is adjusted to 
the change of the electronic states during the transition. 
In the opposite limit the adiabatic approximation is valid 
and FC processes dominate. An estimate for the weight 
of the FC component is 



EX cG 



g,P ^ex,l3' I 



(5) 



which is equal to unity in case of zero mixing and zero 
in case of maximal mixing. The weight of the adiabati- 
cally connected (AC) transition Iac = 1 ^ Ifc is defined 
accordingly. The non diagonal matrix element M is pro- 
portional to the square root of a with a coefficient 77 of 
the order of unity. In the strong coupling regime, as- 
suming that AE = LUi, and np w AE {jip ^ 1), one 
gets 



'FC 



[1 + 4:{Trnp/ncf 



(6) 



where r^p = 1/AE and Tic = l/D. For 77 of the or- 
der of 1, one obtains a robust qualitative description of a 
rather fast switch from AC- to FC-dominated transitions, 
when Ipc and Iac exchange half of their weights in the 
range of a from 7 to 10 (see inset of Fig. 2d). The phys- 
ical reason for such a quick change is the faster growth 
of the energy separation AE ^ o? compared to that of 
the matrix element M ~ a^^^. This switch has nothing 
to do with the self-trapping phenomenon where cross- 
ing and hybridization of the ground state and an excited 
state occurs. The AC-FC switch is a property of transi- 
tions between different states and is related to the choice 
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